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^ ' Abstract: 

Q_i' The process of exclusive heavy vector meson photoproduction, 7p — > ^p, is studied in 
^ ■ the framework of QCD factorization. The mass of the produced meson, ^ = T or J/\l/, 
provides a hard scale for the process. We demonstrate, that in the heavy quark limit and 
at the one-loop order in perturbation theory, the amplitude factorizes in a convolution of a 
^ I perturbatively calculable hard-scattering amplitude with the generalized parton densities and 
' the nonrelativistic QCD matrix element (Oi)y. We evaluate the hard scattering amplitude 
at one-loop order and compare the data with theoretical predictions using an available model 
for generalized parton distributions. 



1 Introduction 



The process of elastic production of heavy quarkonium in photon-proton colhsions, 



was studied in fix target ^ |2] and in HERA colhder experiments both for the case of a real 
photon in the initial state (photoproduction) jHl IH El El E! and for the case when the meson 
is produced by a highly virtual photon (electroproduction) |H1 E! • The primary motivation 
for the strong interest in this process (and in the similar process of light vector meson 
electroproduction) is that it can potentially serve to constrain the gluon density in a proton. 
On the theoretical side, the large mass of the heavy quarks provides a hard scale for the 
process which justifies the application of QCD factorization methods that allow to separate 
the contributions to the amplitude coming from different scales. 

The first step in this direction was made by M. Ryskin ^U] who expressed the amplitude 
of exclusive heavy meson production in terms of the gluon density and, accordingly, predicted 
that the cross section, which is proportional to the square of the gluon density, grows fastly 
with energy. Electroproduction of light vector mesons was studied later in where it was 
shown that in this case the amplitude factorizes in terms of a perturbative hard scattering 
coefficient function and nonperturbative quantities: a meson distribution amplitude and a 
gluon density in a proton. Again, an increase of the cross section with energy was predicted. 
The data from HERA appear to be in accord with these predictions. 

The early approaches to factorization in exclusive vector meson production fOl were 
based on the use of leading double ln(l/x) InQ^ approximation and were designed for the 
description of the process at high energies (in the diffractive or small x kinematics). Later 
on it was understood that in the scaling limit, — ^ oo and x = Q'^/W^ fixed, deeply virtual 
meson electroproduction [T2] and Compton scattering (DVCS) [T21II11II3] processes may be 
studied within the QCD collinear factorization method. The proof of factorization for meson 
electroproduction was provided in jTHj. Due to nonvanishing momentum transfer in the t— 
channel the amplitude of this deeply virtual exclusive process factorizes in terms of general- 
ized parton distributions (GPDs) rather than the ordinary parton densities which enter the 
QCD description of inclusive deep inelastic scattering and the other hard inclusive processes. 
GPDs extend the forward parton distributions and the nucleon electromagnetic formfactors 
to the nonforward kinematics of the electroproduction processes, they encode much richer 
information about the dynamics of a nucleon than the conventional parton distributions. 
This additional information can be presented e.g. in terms of spacial distributions of energy, 
spin ... within a nucleon [TTl UHl CH] • By now the studies of deeply virtual exclusive processes 
and GPDs have developed in a very dynamical field, for recent reviews see, [201 

Another QCD approach to exclusive meson production at high energies is related to 
k±— (or high energy) factorization fl^ I23j. it is based on the BFKL method |211l2ni- In 
this scheme large logarithms of energy In (1/x) are resumed and amplitudes are given by an 
overlap integral of the k± dependent (unintegrated) gluon density and the hard scattering 
kernel. High energy factorization can be formulated also in terms of color dipoles [211 EZj • 
Although these approaches to hard diffractive processes are very promising, their firm foun- 
dation, unfortunately, remains limited due to the use of the leading In (1/x) approximation. 



7P — > Vp , where ^ = J/\E' or T , 
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A few years ago the BFKL formalism was extended to the next-to-leading order |2HI; but the 
generalization of the kj_— factorization scheme or the dipole approach to this order remains 
still a matter of debate. An extended overview of different approaches to this problem may 
be found in 

Most of the theoretical studies HaEBEaESlEllESlEniEZlEH! of process (mH) were 
performed in the framework of kj_— factorization or dipole approaches. Despite of great 
progress and evident success in the description of the data the theoretical uncertainties 
remain poorly understood. In particular, it is believed that the account of skewedness, 
i.e. the effect of different parton momentum fractions, is very important for the kinematic 
range available in the experiment. But since this effect is beyond leading In (1/x), its model 
independent implementation into k±— factorization scheme or dipole formalism remains a 
challenge for theory. 

In this paper we study process in the heavy quark limit in the coUinear factorization 
approach. The physics behind collinear factorization is the separation of scales. The mass of 
the heavy quark, m, provides a hard scale. A photon fluctuates into the heavy quark pair at 
small transverse distances ~ 1/m, which are much smaller than the ones ~ 1 /A related to any 
nonperturbative hadronic scale A. We will show by explicit calculation that to leading power 
in 1/m counting and one-loop order in perturbation theory the amplitude is given by the 
convolution of the perturbatively calculable hard scattering amplitude and nonperturbative 
quantities. The latter are gluon and quark GPDs and the nonrelativistic QCD (NRQCD) [HH] 
matrix element {Oi)v which parametrizes in our case an essential nonrelativistic dynamics 
of a heavy meson system. This means that two firmly founded QCD approaches, namely 
collinear factorization and NRQCD, can be combined to construct a model free description 
of heavy meson photoproduction which is free of any high energy approximation and may 
be used also in the kinematic domain where the energy of the photon nucleon collision, W, 
is of order of the meson mass, M. We evaluate the hard scattering amplitude at next-to- 
leading order. This allows to reduce the scale dependence, which is especially important 
at high energies, since in this case (i.e. in the small x region) the dependence of the gluon 
distribution on the scale is very strong. 

The factorization theorem ^H] for meson electroproduction expresses the amplitude in a 
form containing a meson light-cone distribution amplitude. Its application to the production 
of a heavy meson is restricted to the region of very large virtualities, ^ m^, where the 
mass of the heavy quark may be completely neglected. In contrast, in photoproduction or 
electroproduction at moderate virtualities the heavy quark mass provides a hard scale and 
the nonrelativistic nature of heavy meson is important. In this case, according to NRQCD 
which provides a systematic nonrelativistic expansion, a factorization formalism must be 
constructed in terms of matrix elements of NRQCD operators. They are characterized by 
their different scaling behavior with respect to v, the typical velocity of the heavy quark. 
In the leading approximation only the matrix element {Oi)v contributes, which describes in 
NRQCD the leptonic meson decay rate [5^ 



Here a is the fine-structure constant and m and e are the pole mass and the electric charge 
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of the heavy quark (ec = 2/3, Cb = —1/3). Equation p.2p includes the one- loop QCD 
correction jlU] and as is the strong coupling constant. 

The leading relativistic correction to the meson decay rate and to the photoproduction 
process p.lj) scales ~ (f^), see jHni- It is expressed through the matrix element of an addi- 
tional NRQCD operator. Since for a nonrelativistic Coulomb system v ~ as, the relativistic 
effect is less important than the one-loop perturbative correction. The relativistic correction 
(~ {v"^)) to the result [TOj for heavy meson production was studied in jlj, see also |32]- 
Despite the fact that ~ 0.2 0.25, the relativistic effect was found to be rather 

small. On the cross section level it amounts to 7% for J/\l/ and it should be even smaller for 
T production. 

We will neglect relativistic corrections and consider the process in leading order 

of the relativistic expansion. In this case all essential information about the quarkonium 
structure is encoded in one NRQCD matrix element. In potential models it can be related 
to the value of the radial wave function at the origin. 



here Nc = 3 for QCD. Due to the relation to potential models this scheme of calculation 
is often called in the literature the static or non-relativistic approximation. However, one 
should notice that using NRQCD it can be improved in a systematic and rigorous way 
calculating relativistic and perturbative corrections. In this paper we will concentrate on 
the one-loop perturbative correction. Our main result is that for this process the coUinear 
factorization method is compatible at one-loop level with the relativistic expansion. This 
allows us to obtain unambiguous predictions. We found that QCD corrections are large. 
They change not only the overall normalization but may affect, also, the predictions for the 
dependence of the cross section on energy. 

Our presentation is organized as follows. In Section 2 we introduce the notations, dis- 
cuss the factorization procedure and give the predictions for the amplitude in leading order 
(LO). Section 3 is devoted to the detailed derivation of the hard-scattering amplitude at 
next-to-leading order (NLO). Our method is similar to the one we used recently for the 
calculation of light vector meson electroproduction in NLO. It is based on the use of disper- 
sion relations and the low energy theorem for the radiation of a soft gluon, the non-abelian 
generalization of the theorem known in QED In Section 4 we present a numerical 

analysis. In the concluding section we summarize and discuss open questions. 

2 Factorization and the amplitude at LO 

The kinematics of heavy vector meson photoproduction is shown in Fig. 1. The momenta 
of the incoming photon, incoming nucleon, outgoing nucleon and the produced meson are q, 
p, p' and K, respectively. In the leading order of the relativistic expansion the meson mass 
can be taken as twice the heavy quark pole mass, K'^ = and M = 2m. The photon and 
nucleon are on the mass shell, = 0, p^ = p'"^ = where is the proton mass. The 
photon polarization is described by the vector e^, (e^g) = 0. The invariant cm. energy is 
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s-yp = {q + pY = We define 

{q-Ar = K' = M\ C=^. (2.4) 

In our case the variable ( has a similar meaning as xbj in the electroproduction process. 
We introduce two light-cone vectors 

= = , n+n_ = 1 . (2.5) 

Any vector a is decomposed as 



a' 



^ — a^n^ + a n'l + a_L , a? = 2a^a — . (2-6) 



We choose the light cone vectors in a similar way as in Ji's notation, namely 



^ 2(1 + 

We are interested in the kinematic region where the invariant transfered momentum, 

t = = - (y^< + ^A^) , (2.8) 

is small, much smaller than m. In the scaling limit the variable ^ which parametrizes the 
plus component of the momentum transfer equals ^ = C/(2 — C)- 

The amplitude of quarkonium bound state production can be derived from the matrix 
element which describes the production of the on-shell heavy quark pair, = ^2 ~ ^"^i 
Q'l + i?2 = with a small relative momentum. The explicit equations providing the projec- 
tion onto quarkonium states with different quantum numbers may be found in [13]. For the 
S*— wave, spin-triplet case, which we are interested in, the procedure corresponds to neglect- 
ing the relative momentum of the pair, qi = q2 = and the replacement of the quark 
spinors by 

1 1^ 

^.fe)u,(gi)-^(^) >v(^ + M). (2.9) 
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Figure 1: Kinematics of heavy vector meson photoproduction. 



Here the indices i, j parametrize the color state of the pair, and the vector ey describes the 
polarization of the produced vector meson, (eyCy) = —1 and {Key) = 0. 

CoUinear factorization states that to leading twist accuracy, i.e. neglecting the contribu- 
tions which are suppressed by powers of 1/m, the amplitude can be calculated in the form 
suggested by Fig. 1: 



Here ^^(xi, /i^^) is the gluon or quark GPD in Radyushkin's notation ^2j; xi and X2 = Xi—C, 
are the plus momentum fractions of the emitted and the absorbed partons, respectively. 

is the hard-scattering amplitude and is the (coUinear) factorization scale. By 
definition, GPDs only involve small transverse momenta, k± < fip, and the hard-scattering 
amplitude is calculated neglecting the parton transverse momenta. Since quarkonium con- 
sists of heavy quarks, it can by produced in LO only by gluon exchange. The Feynman 
diagrams which describe the LO gluon hard-scattering amplitude are shown in Fig. 2. The 
contribution of the light quark exchange to quarkonium photoproduction starts in coUinear 
factorization at NLO, it is shown in Fig. 3. Since in this paper we consider the leading 
helicity non-flip amplitude, in eq. ()2.10p the hard-scattering amplitudes A^jf{xi,fip) do not 
depend on t. The account of this dependence would lead to the power suppressed, ^ t/m, 
contribution. 

The momentum fraction Xi, < Xi < 1, is defined with respect to the momentum of 
the incoming proton. It is convenient to introduce the variable x, —1 < x < 1, which 
parametrizes parton momenta with respect to the symmetric momentum P = [p + p')/^. 
The relation between the different variables is 



M 




p=9,g,g 



^ / dxiA^j^{xi,fip)J^^{xi,t,fip) 



(2.10) 



(2.11) 



Xi = 



X2 = 
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Figure 2: The hard-scattering amplitude at LO. 



In terms of symmetric variable x the factorization formula reads 

1 



M 



Airy/Ana eq{eye^) f (Oi) 



1/2 



J dx [ Tg{x, F^{x, e, t) + T,{x, oF^'^ix, e, t) ] 



(2.12) 



q=u,d,s 



Here the dependence of the GPDs and the hard-scattering amplitudes on /i^ is suppressed 
for shortness. In the quark contribution the sum runs over all light flavors, see Fig. 3. 

GPDs are defined as the matrix element of the renormalized light-cone quark and gluon 
operators: 



z=An_ 



2(Pn„) 

F'{x,U) 
1 

~ 2(Pn_) 



W^ix, e, t) u{p') i.u{p) + £\x, t) u{p') 



2tx 



2 



(2.13) 



2mjv 



(2.14) 



XlP+ 



X2P+ 



P 




P 



Figure 3: The light quark contribution to heavy meson photoproduction. 



In both cases the insertion of the path-ordered gauge factor between the field operators is 
imphed. In the l.h.s. of eqs. (I2.13|) . (j2.14j) the dependence of GPDs on the normalization 
point /ii? is suppressed for shortness. In the forward limit, p' = p, the contributions propor- 
tional to the functions S'^{x,^,t) and S^{x,^,t) vanish, and the distributions 7P{x,^,t) and 
T-C^{x,C,,t) reduce to the ordinary quark and gluon densities: 



Note that the gluon GPD is an even function of x, 7i^{x,^,t) = H^{—x,C,,t). 

The definition of the gluon distribution ()2.14|) involves a field strength tensor and, there- 
fore, is valid in any gauge. But to evaluate the gluon hard-scattering amplitude, it is con- 
venient to consider the light-cone gauge n_A = 0. In this gauge the parton picture which is 
behind the collinear factorization formalism appears at the level of the individual diagram. 
One can calculate the contributions of each gluon diagram separately by considering photon 
scattering of on-shell gluons with zero transverse momentum and the physical, transverse, 
polarizations. These gluonic amplitudes have to be multiplied by the light-cone matrix 
element of two gauge field operators, which has the form [T^ 



Here a, b are the gluon color indices, (7^^, = g^^^, — nj^^n_i, — n_fj_n^^. The ie prescription for 
the poles in the r.h.s. of eq. (j2.16j) is important since corresponding singularities lie within 



W{x,0,0) 
H'?(a;,0,0) 
W{x,0,0) 



q{x) for X > , 
— g(— x) for X < ; 
X g{x) for X > . 



(2.15) 
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the integration domain and contribute to the imaginary part of the amphtude. In simple 
terms the sign of ie can be understood in this due to the substitution s s + ie, 

or ^ — >■ ^ — ze. But one should notice that such an argumentation may not work for more 
complicated processes which have in their physical regions the absorptive parts in variables 
other than the energy. For an example and an extended discussion of this issue see [IH]. In 
the case of meson photo- and electroproduction the correct sign of ie is given by eq. ()2.16|) . 

The gluon and the quark hard-scattering amplitudes Tg{x,^) and Tq{x,^) describe the 
partonic subprocesses 

= ^-/G~*{QQ)G (2-17) 

and 

A = A9^{QQ)g' (2-18) 
respectively. Here Q and q denote the heavy and light quark. 

rp ( .X ^ . f x-^ + ie 

T,(x,0 = a(^^^) • (2.19) 

In the first relation the factor ^/{{x — C,+ie){x+C,—ie)) in front of the gluon amplitude comes 
from the parametrization of the gluon matrix element in the light-cone gauge eq. ()2.1(i|) . 

Partonic amplitudes depend on two independent dimensionful variables, the partonic 
subenergy s = Xis and the meson mass = (s. Being dimentionless quantities the 
partonic amplitudes can be expressed as a function of the ratio 

y=i^ = ^ = ^. (2.20) 

This convention is adopted in eq. (j2.19|) . 

Another Mandelstam variable for partonic subprocess is m = — s = —Xis. The 
exchange between the two channels, s m, corresponds to the replacements xi ^ —X2, or 
y ■ir^ —{1 + y) , or X ^ —x. Hard scattering amplitudes and GPDs possess definite symmetry 
properties which are closely related to charge conjugation invariance. A photon and a vector 
meson have the same C— parities, which selects C— even exchange in the t— channel. For 
gluons only a C— even GPD exists at leading twist, which is an even function of x, as 
thus also the gluon hard-scattering amplitude is even in x, Tg{x,C,) = Tg{—x,^). For the 
quark there exist both C— even and C— odd GPDs, and F'^ has no definite symmetry under 
the exchange x ^ —x. But since the quantum numbers of the photon and vector meson 
select the C— even exchange in the t— channel, the quark hard-scattering amplitude obeys 
Tq{x,^) = —Tg{—x,^). Therefore only the C— even (singlet) component of the quark GPD, 
pqi+) ^ F''{x,^,t) - F«(-a;,C,t), contributes to (ET^ . 

Next, we have to evaluate the partonic amplitudes Ag and Aq. We will use the di- 
mensional regularization method, with D = 4 + 2e dimensions, in order to regularize the 
ultraviolet (UV) and infrared (IR) singularities which appear at the intermediate steps of 
the calculation. 
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At lowest order there exists only the gluon contribution. Ag is given by 6 tree diagrams 
shown in Fig. 2. A simple calculation gives the result 

4°)(y)=«5(l + e), (2.21) 
Af\y)=0. (2.22) 

The factor (1 + e) in eq. (j2.2H) appears since in D dimensions there are 2(1 + e) transverse 
polarizations of a gluon. It has to be substituted by 1 when A^^^ is inserted in eq. fl2.19j) 
to get the LO result for the gluon hard-scattering amplitude. However, we kept this factor 
in eq. ()2.21|) since it is important for the correct subtraction of the collinear and ultraviolet 
counterterms in the NLO amplitudes. 



3 The hard-scattering amplitudes at NLO 

At LO the gluonic amplitude is a constant, it is a tree amplitude which has no singularities. 
At NLO the one-loop gluon and quark partonic amplitudes develop a branch cut singularities 
along the lines [0, -|-oo) and (— oo, —1] in the complex plane of variable y, see Fig. 4. We 
will use a method based on the dispersion representation in order to simplify the calculation 
of these one-loop amplitudes. Deforming the integration contour as shown in Fig. 4 one 
arrives at a representation of the amplitude which allows to reconstruct it as a function of 
the variable y from its discontinuities along the cuts [0, +oo) and (— oo, — 1]. Thanks to 
the symmetry properties of the partonic amplitudes discussed above the contribution of the 
brunch cut (— oo, —1] to the dispersion integral may be expressed in terms of the discontinuity 
at [0, +oo). 

We will start with the quark contribution, then we present the more complicated calcula- 
tion of the gluonic amplitude. After that we discuss the renormalization and the subtraction 
of the collinear singularities which lead, finally, to the finite results for the hard-scattering 
amplitudes at NLO. 

3.1 The quark contribution 

The dispersion representation for the quark NLO amplitude A^c^\y) reads 

oo 

A^\y) = i / d./™^<"(.) - , (3,23) 



Here ImA^q\z) stands for the imaginary part of the quark amplitude in the s— channel 
of the quark subprocess. Using the crossing symmetry property, Aq{y) = —Ag{—1 — y), 
the contribution of the n— channel discontinuity was expressed in terms of I'mA^q\z), it 
is given by the second term on the r.h.s. of eq. ()3.23|) . For the quark amplitude one can 
use the unsubtracted dispersion relation, eq. ()3.23j) . A''q\z) ~ const at large z, but due 
to cancellation between s— and u— channel contributions the sum of two terms in the 
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Figure 4: The analytical properties of the partonic amplitudes at NLO in the complex plane of 
y = X2K- 





Figure 5: The s- channel cut diagrams for the quark amplitude. 



brackets vanishes at large 2; as ~ 1/2;^ while each individual term vanishes as 1/ z. Thus the 
dispersion integral is convergent at the upper limit. In other words a subtraction constant 
is not compatible with the symmetry properties of the quark amplitude. 

Among the 6 diagrams which contribute to the NLO quark amplitude only 4 diagrams 
have a discontinuity in the s— channel. They are shown in Fig. 5. It is sufficient to calculate 
the first two diagrams which contain a cut of the light quark and the heavy antiquark lines. 
The line of the heavy quark in these diagrams is not cut since it enters directly into the 
meson vertex and, therefore, is effectively on the mass shell. The other two diagrams in 
Fig. 5 describe the heavy quark cut, their contribution is identical to that one of the first 
two diagrams. 

We present the quark amplitude in the form 



alCp 



(47r)i+T(l + e) V /i 



Am/ 



(3.24) 



here r(. . . ) is the Euler gamma function and Cp = (N^ - l)/(2iVc) = 4/3 is the color factor, 
/i is a scale introduced by dimensional regularization. Calculating the imaginary part we 
find 



TT 



Imlqiy) = 2 



l + 2y 



2yl 



1 



1 + y e l + 2y 



3 + 8y{l + y) 
Ayil + y) 



ln(l + 2y) 



(3.25) 
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Then, inserting this equation into the dispersion integral ()3.23|) we obtain the following 
expression for the quark amplitude. 



' '\l + y y J 24y{l+y) l + 2y 

Jn{-y) + ln(l + y) , ^ fln\-y) ln\l + y) 



where 



i + 2y ' ^ V 1 + 1/ y 

3-4y + 16y{l + y) ^. ^ ^ , 7 + % + l%(l + y) 

2RrT^^ ^^2(1 + 2y) ^-^^^-^ L^,i-l-2y), (3.26) 



L^2(2) = -y" yln(l-t). (3.27) 



3.2 The gluon contribution 

The analysis of the gluon contribution follows the same lines as for the quark case. However, 
one has to take into account that the gluonic amplitude is symmetric under crossing, Ag{y) = 
Ag{—1 — y), and that the asymptotics of A^g\y) at large y is A^g\y) ~ y. Therefore we 
need a dispersion representation of A^t^'^ (y) with one subtraction. It is convenient to perform 
this subtraction at ?/ = 0, the point where the second gluon carries zero energy, since the 
calculation of the amplitude in this point may be considerably simplified making use of a 
low energy theorem for the radiation of a soft gluon. The dispersion representation for the 
gluonic amplitude reads 

oo 

Af\y)-A^^\0) = - [ dzImA^^\z) - , ^ • (3.28) 

^ 3 Ti J ^ \z[z-y) [z + y)[z + y + l)J 



The second term in the brackets represents the contribution of the u— channel cut. Due 
to cancelation between the s— and the u— channel contributions the term in the brackets 
vanishes as ~ l/z^ rather than as ~ l/z^ which makes the dispersion integral convergent. 
Therefore one needs only one subtraction, not two. This can also be expressed in the 
following manner: the term linear in y of the subtraction polynomial is absent, because it is 
not compatible with the symmetry property of the gluonic amplitude. 
It is convenient to introduce the auxiliary quantity Tgiij) defined by 

4"M- ..-..± , i'^SUv)- (3.29) 



(47r)i+T(l + e) V /"^ 

The imaginary part of the gluonic amplitude may be represented as sum of three different 
contributions 

Imlg{z) = ImXf^\z) + ImXf'^\z) + Imlf^^z) . (3.30) 
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Xip+ 



Figure 6: The contribution of the QQ intermediate state to the gluonic amphtude. 

Here ImX'p^\z) represents the sum of 10 diagrams having a QQ cut in the intermediate 
state, see Figs. 6 and 7. Imlg^^\z) gives the contribution of the 24 heavy quark gluon cut 
diagrams and Imlg^^\z) is the contribution to the imaginary part coming from the 24 cut 
diagrams with the heavy antiquark and the gluon in the intermediate state shown in Figs. 8 
and 9. The latter two contributions are equal, 

ImXf3\z) = Imlf3\z) , (3.31) 

therefore it is enough to calculate only one of them, say, ImXf^^\z). We define two contri- 
butions 

UV) - MO) = + 2Xf , (3.32) 



in accordance with eq. ()3.28|) . the decomposition of the imaginary part eq. ()3.30|1 . and 
eq. (HSH). 

3.2.1 QQ- and Qg-cut contributions 

The calculation of the QQ-cnt diagrams shown in Figs. 6 and 7 gives 

-JmjfQ(y) = (y)^ef (y), (3.33) 



where 



arctanh./-^ / q q \ 



Here for shortness we denote two independent color structures by 



C 1 

Ci = Cf , C2 = Cf - ^ = . (3.35) 
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Inserting this result to dispersion integral ()3.28p we obtain 



QQ 3 + 2i/(l + y) ^J-y{l+y) f7 



' 3-%(l + y) 3 + y{l + y){9-y{l + y)) 



+ ^ 1 + % 3 + % 

zco-'^ — ; : — arctan i / 1 arctan 

yil + y) \l + y \l + y y 



(7 + 4y)ci - 2 — C2 



2y(l + y) V 1 + 2/ 



arctan^ \ / q _l r., _l o.,2 



(3-4j,)c,-2i±^'^c.l . (3.36) 



21/(1 + !/) 

Some words about the calculation of integral ()3.28p for X^*^ are in order. Since ImX^'^^z) ^ 
z'^~^l'^ at small 2, the contribution of the region 2; < 5 (where 5 ^ 1) to dispersion integral 
is of the order 



/ 



_3 5^-2 , 2 

7f 



dzz'—^= 7^. (3.37) 





However, this contribution to X^'^, which is singular for 5 — > cancels with the one coming 
from the region 2 > 5 and we arrive at the finite result given by eq. (j3.36p . 

The appearance of integrals like (|3.37p is related to a phenomenon well known in quarko- 
nium physics. The gluon exchange between the nonrelativistic quark pair contains the 
Coulomb like instantaneous contribution. In the NRQCD formalism its contribution has 
to be subtracted from the hard part of the amplitude. Let us discuss the corresponding 
counterterm. 

In a frame where the QQ system is at rest the momenta of the heavy quarks are: 

gi = (m + £,p), q2 = {m + E - e,-p) , (3.38) 
where E denotes the nonrelativistic energy of the pair. The LO amplitude has the form 

M^^ = C [ dp^ip) . (3.39) 



Here C is some factor and "^{p) is the nonrelativistic wave function of the QQ system in 
momentum representation. The integral ()3.39j) is proportional to the value of the wave 
function at the origin 

rfp^(p) ~ Rs{0) . (3.40) 
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Now consider the as correction. The momenta of the quarks after the gluon exchange 

are 

q[ = {m + e\'p'), q'^ = (m + E - e' , -p') . (3.41) 

For the nonrelativistic system the energy and the momentum variables scale as: E, e, e' ~ 
mv"^; IpI, \p'\ ~ mv. With NLO accuracy the amplitude can therefore be written as follows 

^NLO [ dp^(p) 1 1 - "J^"", [ dp' — + O(asv^) 1 . (3.42) 

The first term on the r.h.s. of ()3.42j) is the LO contribution, the second and the third terms 
represent the NLO correction. The later is finite at w — > 0. The second term of ()3.42j) scales 
~ asCp/v, it comes from the instantaneous Coulomb exchange. Eq. ()3.42|) can be easily 
derived considering the integral over the loop momentum, d^+^^g^ = de'dp', and using the 
nonrelativistic limit for the quark propagators. After integration over the loop energy e' we 
arrive at the expression given above for the Coulomb contribution. It can be recognized as 
the exchange potential responsible for the formation of a nonrelativistic meson bound state. 
Indeed, the Schrodinger equation in momentum representation reads 

(E - '-Mp') = - "f' I dpjM^, . (3.43) 
^ m' ' 27r2(27r)2^ 7 {p-p'Y ^ ' 

The second term in Eq. ()3.42|) . the Coulomb counterterm, integrated over dp produces the LO 
contribution, C J dp' '^{p'), which is already taken into account in the first term. Therefore 
the Coulomb counterterm has to be subtracted from Eq. ()3.42|1 . After that one can put the 
quark pair on the mass shell; v^E 0, and p — > 0. 

The advantage of using dimensional regularization is that the quark pair may be put on 
the mass shell even before the subtraction of the Coulomb counterterm. Since at = 0, and 
p = the Coulomb counterterm becomes the scaleless integral, ~ / dp' /p'^ 0, it has to 
be put equal to zero according to the rules of the dimensional regularization method. That 
means that the Coulomb counterterm is zero in this scheme.* The price to be paid for the 
simplification is the appearance of integrals like (j3.37p . They have to be treated as described 
above. We encountered integrals of this kind also in the calculation of 2rg(0). 

Now we proceed to the calculation of 2^^. The imaginary part related with the Qg-cvX 
presented in Figs. 8 and 9 reads 

i/mjfny)=(Y^)^ef(,), (3.44) 
*We are grateful to Maxim Kotsky for the discussion of this issue. 
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Figure 7: The left and the right effective vertices for the QQ-cnt. 



where 



l + 2y{l + y) fci-C2 



1+y 
5ci ci 



e y 4 

'3ci - 4c2 



5Ci - 4C2 . , s 3C2 , 3Ci - C2 

4 ci - C2 y + — + 

y 2 1+?/ 



4(1 + 2?/) 4(l + 2?/)2 V 2 

Ci 3C2 Ci - 2C2 



+ A{ci-C2)y 



9ci - 22c2 5ci - 2c2 



2(1 + 2?/) 4?/2 4(1 + ?/) 



2 , ln(l + 2?/) - -(ci - C2)(45 - 27r2)e . (3.45) 



Expanding Q^^{y) in e one needs to keep, in the hmit of small ?/, the terms which are up 
to linear in e, since in the dispersion integral ()3.28j) they produce the contribution ~ e°. 

Calculating the dispersion integral with ImXg'^^^ we obtain 



+ + 8(1 + 2?/(i + ?/))(-^"^-^) ^"^^ + y^- 



4e 



Cl 

4 



C2 



1 + 2/ 
3 + 7?/(l + ?/) 
22/(1 + 2/) 



2/ 



TT 



2 + ?/(l + ?/)(43 + 100?/(l + ?/)) 8 + ?/(l + ?/)(47 + 61?/(l + ?/)) 

Cl o/. . C2 



Cl- 



96?/2(l + y)2 

l + 2?/(l + ?/)(5 + 14?/(l + ?/)) 



+ C2 



48?/2(l+?/)2 
l + 2?/(l + ?/) 



2?/(l + ?/)(l + 2?/)2 2?/(l + ?/) J 

'\n\-y) \n\l + y) 



ln(2) 



+ 2(ci-C2)(l + 2?/(l + ?/)) 

V 1 +2/ 2/ 
+ ai(?/) ln(-?/) + ai(-l - ?/) ln(l + y) 

+ a2(2/)i:i2(l + 2?/) + a2(-l - y)Li2{-l 



22/), (3.46) 
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Figure 8: The contribution of the Qg intermediate state to the gluonic amplitude, 
where the functions ai and a2 are given by the following expressions: 



Eqs. ()3.36l3.46p define the r.h.s. of eq. ()3.32|) . To finish our consideration of the gluon 
contribution one still needs to evaluate Xg(0), the one- loop amplitude describing the emission 
of a soft gluon. 

3.2.2 The emission of a soft gluon 

The idea of our method is inspired by the famous result of Low known as low energy 
theorem for radiation of a photon. The arguments of Low may be used to constrain the 
amplitude describing the emission of a soft gluon. In a non-abelian case, due to the confine- 
ment phenomenon, the corresponding result has not such a fundamental meaning as in QED. 
Nevertheless, it can be useful, for problems treatable by perturbative methods. We will first 
explain the essential steps of our approach for a simple example, namely the calculation of 
the LO gluonic amplitude ()2.21|) . Then we proceed to the evaluation of Xg(0). 
Let us consider the gluonic process 







l{q)G{x^p)^V{K)G{x2v) 



(3.49) 
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Rl i?2 R3 

Figure 9: The left and the right effective vertices for the Q^-cut. 

at LO in the limit when the emitted gluon is soft; X2 — > 0, xi — > (. With respect to the soft 
gluon the diagrams in Fig. |21may be divided into three groups. In diagrams a) and b) the 
soft gluon is radiated from the on-shell quark line, in diagrams c) and d) it is attached to 
the on-shell antiquark line, whereas in diagrams e) and f) the soft gluon is emitted from the 
virtual antiquark and the virtual quark lines, respectively. In the first two cases the quark 
propagator attached to the soft gluon vertex is close to the mass shell. We call them pole 
contributions, contrary to the third non-pole case which describes the emission of the soft 
gluon from the internal part of the process.^ Our idea is to calculate the amplitude of the 
process (I3.49p in the soft gluon limit considering the pole contributions only. Below we will 
show how using gauge invariance the non-pole contributions may be derived from the pole 
ones. 

Neglecting the proton mass and one has 

W 

p=il + OWn+, g = ___n_, K = q + Cp. (3.50) 

For the photon polarization vector we choose the gauge (e^p) = 0, hence = e^. Since we 
are interested in the helicity non-flip amplitude the meson polarization vector can also be 
chosen transverse, ey = Cy. 

^Duc to color neutrality of the two gluons in the process (|3.49() the emission of gluon G{x2p) from the 
on-shell line of gluon G(xip) is forbidden, thus in our case there is no gluon pole contribution. 
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For the process ()3.49|) in coUinear kinematics it happens that for the pole contributions 
a pole factor Xjx^ coming from the denominator of the quark propagator is compensated by 
the factor X2 from the nominator. This means that contributions of both the pole and the 
non-pole diagrams are regular at 0:2 ^ and that both classes of diagrams contribute to the 
amplitude on equal footing. However in order to apply our method we need to have a pole 
factor in the pole contributions. For this purpose we change the kinematics of the process 
fl3.49|) slightly away from the coUinear one introducing the small transverse component to 
the momenta of the photon and the soft gluon 



q' = q + k^, X2P —^k = X2P + 



(3.51) 



Note that this replacement makes the photon and the soft gluon lines slightly virtual, q'^ = 
k"^ = k\. But this effect is quadratic in k±_ and, therefore, it is small and can be safely 
neglected, as we will always do below. 

The change of the photon momentum (j3.51|) leads to the following replacement in the 
expression for the photon polarization vector 



S TTI^P: 



{pq) 



ieW) 



0. 



(3.52) 



We denote the polarization vectors of the gluons with momenta xip and k by el and el 



(e» = 0, {elk) = 0.^ 



(3.53) 



and choose a gauge such that (e^ q) = (e^ g) = . Thus, the polarization vector of the first 
gluon is transverse, = e^^, whereas the polarization vector of the soft gluon contains 
both a transverse and a longitudinal component, is transverse only in the coUinear limit: 



e'^ at k I 







Let us consider one of the pole diagrams, say, diagram b ). Its contribution to the gluonic 
amplitude reads 



A, 



DSp 



/{/2+ 1^ - xifi+m 



^» {K/2 + kf - m2 ^3 + k- xipf 



(3.54) 



here D is some factor which is irrelevant for our argumentation. The first propagator on 
the r.h.s. of (j3.54j) is the propagator of the quark attached to the soft gluon vertex. Its 
denominator, {K/2 + k)"^ — = (kK), vanishes in the soft gluon limit. In accordance with 
the nominator of this propagator we define two contributions 



A, 



A add 



A add I AZ 
^b) "T ^b) 1 



= DSp 
DSp 



where 
^72 + m 



'{K/2 + kf 



^/2+'l^-Xi^ + m 

rn^ ^ {K /2 + k — xipY — m 

. f[/2+^- xi]^ + m 



(3.55) 



_3 {K/2 + fc)2 - m? '^^ {K/2 + k- XipY 



, fy{^ + M) 

^'fy{^+M) 
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Commuting in the first term the factors ^ and [^f^jl + vti) we obtain 



{kK) 



3{K/2 + k-xipf -m? 



(3.56) 



The trace on the r.h.s. of eq. ()3.56|) vanishes linearly in k. Indeed, using the properties of 
the polarization vectors discussed above it is easy to see that 

Similarly, for the second term we obtain 

= -^^^ Sp [f^ ^ 4\ < fy ^] + 0{k) . (3.58) 

A'^'^ vanishes in the coUinear kinematics since at A;^ ^ 0: k ^ X2P, — e^"*" and [GgK) . 
However, for k^ both contributions to Ai,) are of the same order. Note that A^^ is finite 
for k± = and X2 — * 0. 

The consideration of the other pole diagrams follows the same lines. We calculated, 
similar to eqs. ()3.57| I3.58j) . the corresponding contributions to A"'^'^ and A^ of each pole 
diagram. The total gluonic amplitude is 

A = ej' "A^ = el' " [A'^f'^ + + A^-^"'*^] , (3.59) 

where the first two terms represent the pole contributions, and the third term stands for the 
contribution of the non-pole diagrams. The later can be obtained from A'jf'^ using gauge 
invariance. Due to current conservation we have 

PA^ = . (3.60) 

Since 

= (3.61) 
by construction, see eqs. ()3.551 I3.58|) . we obtain 

In its turn Aff'^ has the form 

<'=(^(fc^), P' = Y.Pt- (3.63) 

The vector P receives contributions from the pole diagrams enumerated by the index i. For 
instance, according to eq. ()3.57|) . the contribution of diagram b) is 

p^) = -^spyir Afv] ■ (3-64) 

19 



Similarly, we denote the contributions of separate pole diagrams to as Af, 

A^ = Y.Al. (3.65) 

i 

From eqs. ()H.fi2j) and (jH.fiHj) we deduce that 

^n-pole ^ _ (3_gg) 

Thus we have shown how in the soft gluon limit the contribution of the non-pole diagrams 
can be derived without explicit calculations. 

Returning to the coUinear kinematics, we have 

A\,^^, = el^^^[Al-P,] |,^_o , (3.67) 

here we used that A"''^'^ vanishes at /c_l — > 0. Thus, the first term in eq. ()3.67|) represents 
the contribution of the pole diagrams in the coUinear limit whereas the second term, ~ P^, 
restores the non-pole contribution to the gluonic amplitude in this limit. 

Finally, to obtain the gluon hard scattering amplitude j^g\y = 0) one needs to perform 
the summation over 2 -|- 2e transverse polarizations of the gluons {,^\\^ = = 
1, . . . , 2-|-2e) in the amplitude of the gluonic process A|fc^^O; and then take the limit X2 — ^ 0. 

Proceeding separately for each pole diagram with these steps, including the summation 
over the gluon polarizations, we find 

Af{y = Q) = Y,D^. D, = Dt + Df\ (3.68) 

i 

where Di stands for the contribution to the gluon hard scattering amplitude of the individual 
diagram. and D'}'^ corresponds to the contribution of the pole diagram i to and 
respectively. After a simple calculation we find 

Dl^ = Dl^ = Dtf' = D^,f = ^^Ii±il , = Dl^ = D:f = = . (3.69) 

Thus we confirm eq. fl2.21|) . Calculating contributions of the non-pole diagrams e) and f) 
directly one can check that the non-pole contribution is indeed correctly restored by jjadd^ 
Although it is very simple this calculation contains all essential points of our method. 

Now we proceed with this method to the evaluation of Tg{0). The one- loop diagrams 
describing the radiation of a soft gluon from the on-shell antiquark line are shown in Fig. E3 
A similar set of diagrams can be drawn for the radiation of the soft gluon from the on-shell 
quark line. Since these two sets of diagrams transform into one another under the charge 
parity transformation it is enough to calculate one of them, say, those in Fig. ^1 and then 
to double the result. 

The results of our calculation of the contributions of individual antiquark pole diagrams 
are summarized in Tables 1 and 2. Using our procedure we obtained for each diagram 
Di, . . . , Du two quantities and Df'''^} 

*The diagrams D2, include the instantaneous Coulomb exchange which we treated in dimensional 
regularization as discussed above. 
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Figure 11: Mass counterterm diagrams which have an antiquark pole in the soft gluon hmit. 

Besides soft and colhnear singularities the one-loop gluonic amplitude contains also ul- 
traviolet poles which have to be subtracted in the on-shell scheme. The full renormalization 
procedure includes mass counterterm diagrams, the renormalization of the heavy quark field 
and the renormalization of the strong coupling constant. The field and the coupling renor- 
malization will be discussed later, together with the factorization of colhnear singularities. 

Here we will consider the mass counterterm diagrams. This can be done in our method 
by considering only mass counterterm diagrams having an antiquark pole in the soft gluon 
limit. They are shown in Fig. ^2 Thus, similar to Df and Df'^'^, we have in Table 2 
two contributions for the diagrams C2 and C4. Below we show that the sets of diagrams 
-D12, -D14, D16 and Di3, D15, D17 together with the mass counterterm diagrams Ci and C3 
add up to two combinations which are gauge invariant. These sets are the separate gauge 
invariant contributions which can be calculated, similar to D^, directly in the colhnear limit. 

At the one-loop level the mass and quark field renormalization constants are equal [Tf] 



Mass counterterm diagrams are multiplied by 6m/m. Let us consider {Di2+Ci^+Di4^+Diq) 
and {Dis + + Di^ + -D17), which represent the one-loop correction to the soft gluon 

vertex 



multiplied by the LO antiquark pole diagrams Bi and B2 shown in Fig. ^1 After a straight- 




(3.70) 




(3.71) 
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Figure 12: LO antiquark pole diagrams. 



forward calculation we obtain^ 



W2 



W3 



Cl 



C2 



3 + 2e 
e(l + 2e) 



3 + 2e 
'e(l + 2e) 



m 



2e 



2e 



(ci - C2) 



3 + 2e 
'e(l + 2e) 



The sum of these contributions equals 

'Cl - C2 



w 



3ci + 2c2 + 0(e) ^ 
e / m 



1 



1 -e 
efl + 2el 



(3.72) 



(3.73) 



(3.74) 



(3.75) 



We found that for the one-loop correction to the soft gluon vertex the contribution ~ 
cancel. What is left is which means that the part ~ as of (j3.71|) is gauge invariant. 

Inserting it into the LO diagrams Bi and B2 we obtained the results presented in the first 
two lines of Table 2. 

Note that in the abelian case ci, C2 = 1, and, according to ()3.75|) . the correction to the soft 
vertex is finite at e — > 0. It corresponds to the contribution of a fermion anomalous magnetic 
moment, a/(27r), which is in accordance with the general statement of Low's theorem in 
QED. In the non-abelian case this contribution has no such clear physical meaning since it 
is infrared divergent. 

Finally, summing all contributions in Tables 1 and 2 and multiplying the result by a factor 
2 (thus taking into account the quark pole diagrams) we arrive at the following expression 
for the gluonic amplitude at X2 = 



, Cl-C2 7ci-C2 

' e2 2e 



3 37r2 
'2 + - 



101n(2) j-C2 ( 5 



57r^ 



21n(2)j . (3.76) 



^Note that wi is finite for fc — > only if the mass counterterm diagram is included. 
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3.3 NLO results 



Let us discuss the structure of singularities of the parton amphtudes. First of all note 
that, according to ()3.32j) . the double pole terms which are present in eqs. fl3.7(i|l and ()3.4(ij) 
cancel. Thus, the gluonic amplitude as well as the quark one contains only single poles in e. 
This means that soft singularities present in the individual contributions cancel in the final 
expressions for the NLO parton amplitudes. What is left are the single poles in e which as 
we show below represent the ultraviolet and coUinear singularities. 

To demonstrate the validity of factorization one needs to check that the ultraviolet poles 
are removed by the heavy quark field and the strong coupling renormalization, and that the 
coUinear poles are absorbed into the quark and gluon GPDs. For this purpose let us recall 
the structure of the factorization formula 



1 

M dx 



ff (x, + T^(^) (x, O) F'{x, e, t) + f (x, OF^'H^, e, t)] , (3.77) 



-1 



where the tilde indicates that the renormalization and the separation of the coUinear singu- 
larities has yet not been performed. The bare hard-scattering amplitudes are 



here A'^'^ is defined by eqs. (jT^ . (jS321), and the NLO gluonic amplitude A'^^ by eqs. (jT^ . 
fl3.32|l . ()3.3(j|l . ()3.4(j|l and ()3.7(jj) . In ()3.78|1 and in some equations below we suppress for 
shortness the ie prescriptions, they are easily restored by the replacement ^ ^ — i£- 

The factorization of the coUinear singularities corresponds to the substitution, in ac- 
cordance with the definition of GPDs, of the bare quantities F'^'^{x,^,t), F^{x,^,t) by the 
renormalized ones. In the modified minimal-subtraction (MS) scheme one has at the one-loop 
level 

F-^{x, t) = F-^ix, e, t, - ^ (i + In ) X 

1 

X j dv [Vqq{x, v)F'''^{v, e, t, /i^) + Vggix, v)F\v, e, t, /i^)] , (3.79) 
-1 

27r \e \fi' 



F^{x,U) = - ^ ( ^ + ln ( ^ ) ) X 



1 



X / dv [Vggix, v)F^{v, e, t, /i^) + Vgg{x, v)F'^'''{v, e, t, /i^)] , (3.80) 



-1 
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where Vqq, Vgg, Vgq, Vqg denote the one-loop evolution kernels. 



+ 7s - ln(47r) 



(3.81) 



•jE is Euler's constant. Inserting ()3.79|) and ()3.8()|1 into eq. ()3.77j) and truncating the series 
at the order a| we found the following coUinear counterterms to the gluon and quark hard- 
scattering amphtudes 



(3.82) 



Af(x,0 



as f 1 



2ti \ e 



In 



/4 



dvfS'Xv^OVgqiv.X 



(3.83) 



Note that, since Tf ^ = for our process, the renormalization of the quark GPD ()3.79|) 
does not generate contributions (~ Vqq, Vqg) to the coUinear counterterms. Calculating the 
integrals ()3.82p . ()3.83|) with these kernels we obtain 



Af(x,0 



at 



T + l + ln(4 



at 



271 \ e 



2i 



NcCg 



X — ^ 



2 



where 



/5o 



lliVc 2n 



f 



3 3 ' 

rij is an effective number of light quark flavors, 

C,(,)^,1.2,)(i^-Ml±l) 



, (3.84) 
(3.85) 

(3.86) 



(3.87) 



For the renormalization of the strong coupling one has to substitute the bare coupling 
constant as by the running coupling asififi) in the MS scheme. 



"5 = «s(At 



R) 



1 + 



asijJ, 



R) 



47r 



f^R 



(3.88) 



This substitution generates the following counterterm to the gluon hard-scattering amplitude 



at 



1 



47r ix-0{x + Ve 
25 



1 + ln 



^2 



/5o 



(3.89) 
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Figure 13: The cross section for T photoproduction; the theoretical predictions at LO for the 
scales \ip = fiji = [1.3, 7] GeV (ranging from bottom to top), and the data are from ZEUS [5] and 
HI IS]. 



To account for the heavy quark field renormalization effect one has to add the counterterm 

Af(x,0=5^2Tf(x,0, (3.90) 

with 6Z2 given in eq. ()3.70|) . 

In the sum of the bare hard-scattering ampHtudes and the counterterms described above 
all poles in e cancel. Thus, we can now take the limit e — > 



and arrive at finite results for the hard-scattering amplitudes: 

' X — C, + ie 



J e^O 



zn 



2i 



(3.91) 



(3.92) 



/,(,) = (ln^-l)(l + 2y) 



In(-y) ln(l + y)\ 2 13(1 + 22/) , 2 In 2 



y 



y 



TT 



+ 



A8y{l + y) l + 2y 



In(-v) + ln(l + v) 



In^(-y) ln^(l + y) 

1 + 2/ y 



3-% +162/(1 + 2/) 7 + 4y + 16y{l + y) ^. , ^ 

-Li2{l + 2y) m~, — ^ Li2{-l-2y), (3.93) 



for the quark, and 



Ay{l + y) 

e 



X — ^ + ie){x + ^ — ie) 



Ay{l + y) 



(3.94) 
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/,(y)=4(ci-C2)(l + 2y(l + y))(^ 



y y i^f 



+ 4(ci-C2) (1 + 22/(1 + 2/)) 



\n\-y) \n\l + y) 

i + y y 



2 + y(l + 2/)(25 + 88y(l + y)) 10 + y(l + y)(7 - 52|/(1 + y)) 

Ci H , 7^5 C2 



Cl 



482/2(1 + |/)2 
1 + 62/(1 + 2/)(l + 22/(1 + 2/)) 



242/2(1+2/)^ 



2/(l + 2/)(l + 22/)2 



+ C2 



:i+22/)^ 



y{^ + y) 



ln(2) 



+ 7r 



A/-y(l + 2/) 
2/(1 + 2/) 



-Cl - 3C2 



2C2 



^/-y{l + y) fl + ^y 
2/(1 + y) 



arctan 



3 + 42/ 

y 



arctan 



arctan^ ^ ^ 
22/(1 + 2/) 



arctan^ y^i^J 
2y(l + 2/) 



.7^, ^ p i + 22/ - 22/2 

(7 + 42/)ci - 2 — C2 

1 + 2/ 

, N , 3 + 62/ + 22/2 

(3 - 42/)ci - 2 C2 

2/ 



+ 2 ai(2/) ln(-2/) + 2 ai(-l - 2/) ln(l + y) 

+ 2 a2(2/)Lz2(l + 22/) + 2 a2(-l - 2/)^^2(-l - 2y) , (3.95) 

for the gluon. ai(2/), 02(2/) are defined in eqs. ()3.47|1 . ()3.48j) . The expressions in ()3.92p - ()3.95j) 
represent the main result of this paper. 

At high energies, ^ M'^, the imaginary part of the amplitude dominates. The 
leading contribution to the NLO correction comes from the integration region ^ <^ |x| <^ 1. 
Simplifying the gluon ()3.95|1 and the quark ()3.93|1 hard-scattering amplitudes in this limit 
we obtain the estimate 



M 



-4i7r2A/47raeg(eye^) ( {Oi)v 



1/2 



X 



X 



«.(a.h)F'K, t) + ^M^^^l^ In ( ^) [ ^F»(., i, t) 



TT 



(3.96) 



Given the behavior of the gluon and the quark GPDs at small x, F^{x,C,,t) ~ const and 
we see from (I3.96p that the relative value of the NLO correction is 
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Figure 14: The cross section of tiie T piiotoproduction; theoretical predictions at NLO for the 
scales hf = fiR = [1.3,7] GeV (ranging from top to bottom), and the data from ZEUS f5| and HI 



parametrically large at small ^, 



Nr 



(3.97) 

The gluon correction in (j3.97|) is negative unless one chooses a value of the factorization 
scale yUi? < m, which is substantially smaller than the kinematic scale M = 2m. The 
quark correction is also parametrically large at high energies. It is expected to be sizable 
since it collects the contributions of all the light quarks and antiquarks. These qualitative 
observations are supported by the numerical analysis. 



4 Numerical analysis 

We assume as values of the quark pole masses: rric = 1.5 GeV, m{, = 4.9 GeV. {Oi)v 
was evaluated using eq. ()1.2p with as = asifin)- For the generalized parton distributions 
we adopt the parametrizations, evolved both in LO and NLO, of ^Hj that are based on 
the CTEQ6 set of forward distributions [IH]. We neglect the contributions proportional to 
S'^{x,^,t) and S^{x,^,t). In the numerical calculations we use LO strong running coupling 
and LO GPDs and NLO coupling and NLO GPDs for LO and NLO observables correspond- 
ingly. 

Let us start with T photoproduction. We calculate with our formulas the forward ampli- 
tude and the forward differential cross section, da/dA^ at = 0. For the dependence 
we assume, in accordance with the measurements at HERA, the simple exponential 

da f da , \ hK2 1 f da , \ , , . 

For the slope parameter we use h = 4.4 GeV~^. 



28 





Figure 15: T photoproduction, NLO prediction for fj.p = fi^ = 4.9 GeV and its decomposition 
into different contributions, see text. 
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In Fig. ^1 the LO predictions for the total cross section of T photoproduction are shown 
as a function of energy, the data points are from ZEUS [S] and HI [H]. The curves correspond 
to different values of the factorization scale fip which is chosen equal to hr. The experimental 
uncertainties are large. We find that for the broad interval of scales, fip = fJ^R = 1.3-i-7 GeV, 
our LO predictions lie within the experimental error bars. The strong dependence of the 
predictions on the factorization scale is related to the well known fact that scaling violation 
is large for small x. At small x the gluon density increases rapidly with growing fip which 
leads to an increase of the LO predictions with fip. In Fig. we present the results 
of the NLO calculations for the same set of scales. For meson production in NLO this 
effect is partially compensated, as it should be, due to the dependence of the gluon and 
the quark hard-scattering NLO amplitudes on ^p, see eqs. ()3.92p - ()3.95p . As a result we 
observe a substantial reduction of the scale ambiguity of the theoretical predictions in NLO 
in comparison with LO. 

The NLO predictions are generally smaller than the LO ones. The reason is twofold. 
First, according to the parametrizations we use, in this kinematic region the gluon GPD 
in NLO is about a factor of two smaller than the gluon GPD in the LO. This is another 
manifestation of the large scaling violation effects at small x. Second, we find, in accordance 
with the estimate ()3.96p . that the part of the gluon NLO hard-scattering amplitude ()3.94p 
~ a| leads to a contribution which is large and has at fip > m the opposite sign as the 
contribution ~ as induced by the Born term of ()3.94|1 . The last statement is illustrated 
in Fig. ^1 where the different contributions to the NLO result for fip = fiR = 4.9 GeV 
are shown. In the left upper panel we present the cross section. Here the curve labeled 
Born represents the results calculated using only the Born term of ()3.94j) . the other curve, 
labeled total, is the cross section calculated with the complete result for the NLO hard- 
scattering amplitudes, including the quark contribution. To avoid misunderstanding, in both 
calculations the NLO GPDs were used. We see that the parts of the NLO hard-scattering 
amplitudes ~ Og make the cross section significantly smaller. On the bottom left and the 
bottom right panels of Fig. we present the decomposition into different contributions 
of the imaginary and the real parts of the NLO amplitude divided by its absolute value, 
Im Ai/\A4\ and Re Ai/\Ai\. The curves labeled total represent the results calculated with 
the complete NLO hard-scattering amplitudes (j3.94j) and (j3.92|) . In these figures the quark 
contribution and the decomposition of the gluon contribution into the Born and the part 
induced by the term~ a| of ()3.94|) are shown separately. The corresponding curves are 
labeled as 'quark', 'Born' and 'gluon'. In the right upper panel we show ImAi/\Ai\ and 
Re A4/\Ai\ together. We see that, despite the fact that the value of a strong coupling 
constant is small at this scale, asifin)/ {^n) ~ 0.033 at = 4.9 GeV, the gluon correction 
constitutes ~ 30% of the Born contribution with the opposite sign. The quark contribution 
is about ~ 15%, with the opposite sign with respect to the Born contribution. Note also 
that at high energies the imaginary part of the amplitude is about twice the real part. 

The reduction of the ambiguity for the theoretical predictions due to a variation oi ^p 
is even more pronounced if one chooses a fixed value of renormalization scale, see Fig. 
In this case the value of the cross section is predicted to be smaller than for equal scales, 
jJ'R = fJ'P, compare Figs. IT^ and IT^ 

To summarize our results for T photoproduction we conclude that the NLO corrections 



30 



b 0-8 
0.6 



0.4 



0.2 



■ ■ " 50 100 150 200 250 

W[GeV] 

Figure 16: The cross section for T photoproduction, and NLO predictions for the scales fJ-r = 
[1.3, 7] GeV, fiR = 5.9 GeV. 

stabilize the theoretical predictions with respect to variation of the factorization scale. The 
NLO corrections are numerically important. They make the NLO cross sections smaller than 
the LO ones, and for the GPD-model used our results seem to lie somewhat below the data. 

For the photoproduction of J/"^ the situation is different than for T production since 
in this case a value of the hard scale, the quark mass, is smaller. The NLO corrections 
are much larger than in the case of T production for the following two reasons. First, the 
value of the QCD running coupling is larger at smaller scales. Second, the value of and, 
consequently, the effective values of x in the factorization formula is about two orders of 
magnitude smaller than for T production. Therefore the effect of the enhancement of the 
NLO correction at small x, see eqs. ()3.96p . ()3.97|) . is much larger for J/\E^. This is illustrated 
in Fig. fTTI where the labeling of the curves is the same as in Fig.^l The data are from E401 

and ZEUS [7j. Note that contrary to Fig.[T3l in the left upper panel of Fig.[T7|the results 
for the forward differential cross section are shown. We see that although the predictions 
for the cross section are in reasonable agreement with the data, the absolute value of the 
NLO correction is very large. Note also that in this case the quark GPD makes a significant 
contribution. The sum of the gluon and the quark NLO corrections is as much as twice 
of the Born contribution and of opposite sign. Therefore in NLO the total amplitude has 
the opposite sign as in LO. Note also that the imaginary part of the NLO amplitude goes 
through zero aX W ^ 25 GeV, which is unnatural. Thus, we conclude that for the J/\l/ 
photoproduction the higher order corrections are not under control. 

5 Summary 

We have shown by an explicit calculation of the partonic one-loop amplitudes that in the 
heavy quark limit the coUinear factorization and the nonrelativistic QCD approach applied 
to quarkonium photoproduction are compatible and lead to the unambiguous predictions 
(|3.92p - (j3.95p for the hard-scattering amplitudes in the NLO. Presumably such a factorization 
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Figure 17: The differential cross section for J/^lJ photoproduction, NLO predictions for fir = fJ-R = 
1.52 GeV, and the data from E401 and ZEUS [7j. The labehng of the curves is the same as in 
Fig.Hni 

scheme can be generalized to all orders of the strong coupling expansion. The study of this 
issue, although very interesting, goes beyond the scope of the present paper. 

The numerical analysis for the T photoproduction shows that in comparison to LO the 
NLO corrections to the hard-scattering amplitudes reduce significantly the ambiguity of the 
predictions related to the choice of factorization scale. The NLO corrections are large, at 
HERA energies they constitute about ~ 40% of the Born contribution at the amplitude level 
and are of the opposite sign compared to the Born contribution. 

Contrary to that, we find that for the photoproduction of J/ip in HERA kinematics the 
magnitude of the NLO correction is about two times larger than the Born contribution. Also 
we observe a very strong dependence of the theoretical predictions on fip. That forces us 
to conclude that at high energies for J/ip photoproduction these corrections are not under 
theoretical control if one works in NLO, i.e. in the collinear factorization scheme truncated 
at the second order of the strong coupling expansion. 

Note that all steps of the dispersion method developed in this paper can be applied 
directly to NLO electroproduction, the process of a heavy vector meson being produced 
by a virtual photon. In this case the calculations may be much more involved due to 
the presence of an additional parameter, namely m/Q. However, for electroproduction the 
photon virtuality shifts the hard scale to the higher values in comparison to photoproduction. 
This gives hope that the factorization approach may be reliable for electroproduction of J/\l/ 
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starting from some, not too high, values of Q. 

We show, see eqs. ()3.96p . ()3.97|) . that convolution of the NLO hard-scattering ampli- 
tudes with GPDs produces contributions which are parametrically enhanced at high ener- 
gies, ~ a|ln(l/^). These contributions originate from the diagrams of partonic subprocess 
with gluon exchange in the t-channel and are related to the s— channel radiation of an inter- 
mediate parton in the wide interval of rapidity, away from the photon fragmentation region. 
In higher orders such radiation generates contributions ~ ^^(a^ ln(l/^))". The k±_- fac- 
torization approach allows to sum this class of logarithmic corrections to all orders in 05, 
consistently with the fixed-order factorization of the collinear singularities [SD]. It would be 
very interesting to perform such studies for heavy vector meson production. We believe that 
a resummation of these contributions will lead to much more stable theoretical predictions 
at high energies. 
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Table 1: Contributions to Tg{0) of diagrams Di, . . . , Du. 



Diagram 




T^add 
-^10 

T~)add 3 
-^10 8 

Df^ 

Dff -(ci-C2)^ 
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Table 2: Contributions to lg{0) of D12, . . . , -D17 and the mass counterterm diagrams. 
Diagram e^^ e° 

(Di2 + Ci^ + Z}l4 + /^16) -i(Cl-C2) Ci(| + ^)-C2(i + if^) 
(Di3 + C3^ + I)l5 + /^17) -(Cl-C2)i 



r^addSm ^ 91n(2) 

'-"2 m 16^1 ^1 8 

Cf— -cil 

raddSm 3 „ r M -U 31n(2) \ 

*-^4 — -16^1 ^lU + ^—j 
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